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SUMMARY 


The flov and heat transfer are analyzed at the reattachment zone of 
two- dimensional separated laminar boundary layers. The fluid is considered 
to be flowing normal to the wall at reattachment. An approximate expres- 
sion is derived for the heat transfer in the reattachment region and a 
calculated value is compared with an experimental measurement. 


INTRODUCTION 


The mechanism of heat transfer through separated regions is very 
complicated and little understood. Chapman in reference 1 first analyzed 
this problem and obtained an estimate of the average heat transfer in this 
region. A survey of the literature reveals that no local heat-transfer 
analysis has been done for a Chapman type separated region. Reference 2 
indicates that the complete physical problem is too complex to be ade- 
quately described by a simple flow model. Existing experimental work 
(ref. 3 ) shows that the maximum heat transfer in a separated region occurs 
at the reattachment point. References 1 and k indicate that one may be 
able to analyze the flow and heat transfer near this point. This paper 
is concerned with heat transfer in the reattachment zone for normal reat- 
tachment of a two-dimensional separated laminar boundary layer. Figures 
1 and 2 show the model considered. 

The flow approaching the reattachment zone is considered inviscid but 
rotational. The viscous effect is assumed to be confined to the boundary 
layer which develops along the x axis . The partition into inviscid and 
viscous regions is justified because, as will be seen subsequently, the 
vorticity in the inviscid region is smaller than that in the boundary layer 
by an order of magnitude. In the following analysis the fluid in the 
inviscid region is assumed to be incompressible and a closed-form solution 
of the flow field is developed. This solution is then used to solve the 
boundary- layer equations for heat transfer along the x axis . 
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NOMENCLATURE 


h total enthalpy 

K function defined by equation (15) 

L length of reattachment zone 

Z length of separated mixing layer 

lw PrZ 

Nu Nusselt number, - r 

(h e - h w )p e 

P pressure 

Pr Prandt 1 number 

q v heat transfer to the vail per unit area per unit time 
R nose radius of a hypersonic blunt body 

P e u e ^ 

Re Reynolds number, 

^e 

s variable defined by equation (13) 

u x component of velocity 

u e streamvise velocity at the outer edge of mixing layer 
v y component of velocity 

Vq v(0, L) 


X 


*0 


distance along the vail from reattachment point 

particular value of x greater than L but still near the 
reattachment point 


Y 

y 

z 


y 

L 

distance in direction normal to vail 

variable defined by equation (13) or function defined by equation 

(A7 ) 


A 


H On v_n 



C^VJl > 


3 


A parameter defined by equations (3) and (23) 

p dynamic viscosity 

i dummy variable 

p density 

cp,X Laplace transforms defined by equations (13) 
stream function 
n vorticity 


1 


Subscripts 


e outer edge of mixing layer 

r average value for reattachment zone (0 < x < L) 
t reattachment point 

v vail 


INVISCID FLOW REGION 


The flcv approaching the reattachment zone shall be considered 
inviscid, incompressible, but rotational. Figure 2 shows the flov model 
studied. For this case, the distribution of the streamlines on the x-y 
plane is given by the equation 

V s * = (1) 


where \|r is the usual stream function, defined as: 


u 



5x 


( 2 ) 


To obtain an expression for ft(\|r), the velocity profile of the 
incoming stream at a distance y = L must be known. The distance L is 
defined such that the incoming stream is unaffected by the existence of 
the reattachment wall for y > L; that is, u = 0 for y > L and u is 
positive for y < L. Actually u will be slightly negative for y > L, 
as the fluid is being entrained into the mixing layer from the separated 
region, and u is positive for y < L because of the wall. The value 
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of L will be obtained a posteriori from the solution of the flow field. 

A study of the mixing layer solution of reference 1 shows that the stream- 
wise velocity distribution can be quite accurately represented by the 
expression 


From the definition of L 


and 


v(x, L) = -Yoe"^ 

u(x, L) = 0 
8u(x, L) 


Sy 


= o 


( 3 ) 

( 4 ) 

( 5 ) 


The vorticity, n(\[r)^ may be evaluated at any boundary for it is 
constant along a streamline. In the present study, the vorticity is 
evaluated along the boundary at y = L. The vorticity distribution cor- 
responding to the velocity profile of equation (3) and to condition (5) 
Is derived in the following manner. From equation (3), 


%=L = f| = Av 0 e-?« (6) 

To express equation (6) in terms of we write equation (3) as 


v 


y=L 


at _ 


dx 


= -v 0 e 


-Ax 


When the above equation is integrated to satisfy the definition of the 
dividing streamline, t(0, L) = 0, (see ref. l' there results the 
relationship , 


A 

5 

6 

l 


t(x, L) = (1 - e-^-) (7) 

Equation (6) is now rewritten with the aid of equation (7) as 

%= L(t) = Av 0 - (8) 

Equation (8) shows the relationship between the vorticity and the stream 
function at y = L. Vorticity, however, is constant along a streamline 
throughout the flow field. Equation (8), therefore, is applicable for 
the entire inviscid region. 
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The inviscid flow equation (l) and the complete boundary conditions 




: < 00 and 0 < y < L as : 


+ V 0 A = 0 

Sy 2 

(9) 

O 

II 

o 

(10) 

o 

ii 

o 

(ID 

L) = ( 1 - e-^x) 

A 

(7) 

Ii(x, y) | < M 

(12) 


where M is an arbitrarily large positive number and equation (12) 
expresses the bounded condition. The solution of the above boundary value 
problem is obtained by the Laplace transform method as follows. 

The Laplace transforms of the stream function are first defined as 


<p( s, y) = Lx{t(x, y)} 
X(s, z) = Ly{cp(s, y)} 


(13) 


The transforms of the boundary conditions (ll), (T)> and (l2) are also 
needed in the course of the solution. These are, respectively, 


cp(s, 
cp(s, L) 

and 

icp(s, y) | < If 


0) = 0 

(lla) 

v o 

s(s + A) 

(7a) 

for all s > 0 

( 12a) 


where If is another arbitrarily large positive number. When the Laplace 
transformation is performed twice in equation ( 9 ) and the boundary con- 
ditions (10) and ( 11a) are applied, the following equation results. 
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x(s, z) = 


A(z) + B( s) - (v 0 A/sz) 


Yr + z‘ 


K 2 = s 2 - A 2 (15) 

and A(z) and B(s) are unknown but particular functions of z and s, 
respectively. An inverse Laplace transformation of equation (l4) with the 
aid of boundary condition (7a) yields 


<*■* ^ - I sln *J 0 [sln K(L - 5)1A(5)dl 


v 0 A sin Ky , T _^ 1 . „ r 

+ — — (l - cos KL) + — sm KC 

sK 2 K 


’ sin K(y - |)]A(|)d| 


voA sin KL /. „ \ 

— (1 - cos Ky) 

sKr 


Now, in the above equation, the integrals which include the unknown 
function A(|) in the integrand will be evaluated with the aid of the 
bounded condition (l2a). Equation (l6) shows that the first term on the 
right-hand side of the equation, l/sin KL, becomes oo when 

K = ± (17) 

i_J 

-where n = 0, 1, 2, . . • • From equations (15) and (l7)> this means 
that l/sin KL becomes 00 when 


The boundary condition (l2a) implies that the function y) must be 

bounded for all values of s > 0. It is, therefore, necessary that the 
function in the braces in equation (l6) be zero when s = N /(nir/L) + 7\ 
in order to satisfy the bounded condition (l£:a) . The integral with the 
unknovn integrand, therefore, becomes 


OaVJI 
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sin ^ (L - |) 


A(|)d| 


njtv 0 f A[l - (~l) n ] 

' (n:n:/L ) 2 nrt/L)^ + A 2 


+ 7( nir/L) 2 + A 2 [^(nrt/L) 2 + A^ + Aj / 

(19) 

Finally, the complete inverse transformation of equation (l6) is obtained 
by finding the residues at all the single poles for s < 0. The solution 
of the boundary -value problem, equation ( 9 ) with its boundary conditions, 
is thus obtained and is 


f(x, Y) 

v 0 L 


co ^ . ( sinh ALY) + 1 " C ^ h — - £ exp(-ALX) 


00 



n=i 


AL(sin httY) 
njt[ (rut) 2 + (AL) 2 ] 


exp 


-7(nrt) 2 + (AL) 2 X 


( 20) 1 


The value of AL is found by satisfying condition (4) with the aid 
of the solution, equation (20) . At y = L, from the continuity equation, 
u, if other than zero, is a monatonically increasing function of x with 
u(0, h) = 0 and a maximum at x = 00. Thus, condition (4) can be satisfied 
for all x if one sets u(°°, L) = 0. Equation (20) can be used to show 
that 


u(°°, L) _ 1 

v Q ~ sinh 9 vL 


( 21 ) 


The velocity components, u(x, y) and v(x, y), were found to remain 
practically constant for |u(°°, L)/v 0 | < 0.01. Thus, condition (4) can 
be considered to be satisfied when u(°°, L)/v 0 — 0.01, and from equation 
( 2 l) AL is found to be 5 . 3 . In subsequent numerical work 

AL = 5*3 ( 22 ) 


is used. 


A study of reference 1 shows that A in equation ( 3 ) can be expressed 
quite accurately by 


A = 


\/ Re 


2.222 1 


( 23 ) 


1 Equation (l) can also be solved by the method of separation of 
variables (see appendix). 
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Thus, for a given Re and 2 , L and A can be readily determined from 
equations (22) and (23). 

The velocity and pressure distributions along the vail, obtained from 
equation (20), are shovn in figure 3* 


BOUNDARY LAYER AND HEAT -TRANSFER ANALYSIS 


To evaluate the heat transfer along the vail by conventional boundary- 
layer theory one must first investigate the effects of the vorticity and 
the enthalpy gradient of the inviscid region on the boundary-layer solu- A 

tion. The absolute magnitude of the vorticity interaction parameter 
(ratio of average vorticity in the inviscid region to that in the boundary 
layer) vas estimated to be less than 0.1 for 0 < X < 1 and Re > 10 4 . 

Classical boundary- layer theory can be applied with good engineering 
accuracy when the interaction parameter is of this order of magnitude 
(ref. 5 )* The effect of an enthalpy gradient in the inviscid region on 
the thermal boundary layer vill be about the same as that of the inviscid 
vorticity on the momentum boundary layer. 

The present study concentrates on the reattachment zones that may 
exist on a hypersonic vehicle. The heat-transfer analysis (an approxima- 
tion) is based on the theory of a highly cooled boundary layer for hyper- 
sonic blunt bodies, developed in reference 6, and the pressure and velocity 
distributions obtained in the preceding section. 

A typical heat-transfer variation near the reattachment point is 
shovn in figure 3* In viev of the drastic variation of the local heat 
transfer vithin the small distance L, vhich is of the order of the mixing 
layer thickness, an average heat transfer in this zone is of greater 
engineering interest than the local heat transfer shovn in figure 3- The 
reattachment zone is defined as that area of the vail along vhich flov 
readjustment takes place and is defined tvo -dimensionally by 0 < X < 1. 

The average heat transfer at the reattachment zone vas calculated for 
several pressures in the separated region and resulted in development of 
the folloving semiemplrical expression 

q. w , r * 2~ 3 /2 Pr' 2 (h t - h v )(o.76 + 1.4ll h) (24) 

for 0.1 < P e /Pt < 0.5- To express the heat transfer in terms of the 
fluid properties at the outer edge of the mixing layer, the folloving 
approximations are made: 


H CTVJI 
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From reference 6 : 




Pt 

Pe 


From reference 1 for Pr close to 1: 


ht 


■L* * 12. = 0.587 


h e - h v u e 

Equation (24-) now becomes, with the aid of equations 

- 1/2 


Nu r SS 


q-v.r ^ 1 
(h e - h w )n e 


/p v / 

0.01*63 Pr 1/3 Re 3/4 [0 


(22) and (23), 

76 + i. 4 n 

(25) 


The pressure ratio E e / P t can be obtained from the assumption of an 
isentropic compression along the dividing streamline of the mixing layer 
as shown in reference 4- . It should be remembered that the area of the 
reattachment zone varies with Re. 


Sometimes, it may be desired to calculate the average heat transfer 
to an area which is a bit larger than the reattachment zones, but includes 
it. From figure 2 and the definition of L one can see that the pressure 
along the vails is essentially constant for x > L. Hence, the average 
heat transfer for L < x < x 0 is readily found to be 


Nu^ 


L 


< X < X Q ) 


0.0926 Pr 1 / 3 Re 3/4 
0.2922 (x 0 /z) 1/2 Re 1/4 + 1 


(26) 


A weighted average of equations (25) and (26) would give the average heat 
transfer near the reattachment point for x Q > L. 


DISCUSSION AND CONCLUDING REMARKS 


There are no experimental data available, to the authors 1 present 
knowledge, which could be used directly to compare the above theory. 
However, a rough comparison can be made with the experimental results of 
reference 3 in which an average heat transfer for x G / l = 0.08 was 
measured for a separated laminar boundary layer reattaching at an angle 
of about 4-5°. The experimental conditions of reference 3 vere used to 
calculate the average heat transfer which was found to be within 10 
percent of the measured value . 
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The present analysis is for the normal reattachment only and therefore 
its comparison with the experimental result for the 45 ° reattachment angle 
lends only rough support. It is doubtful, however, that the variation of 
the reattachment angle between ^5° and 90° would change the heat -transfer 
result drastically. 

It is interesting to investigate the dependence of the average 
laminar heat transfer on Reynolds number. Equations ( 25 ) and ( 2 6) indi- 
cate that the Nusselt number varies with Re 3 ' 4 in the reattachment zone 
and with a slightly smaller power of Re outside this zone. Thus, the 
average heat transfer for a given x D > L varies between Re 1/2 and 
Re 374 , the exact value depending on x Q and Re. It Is worth noting that 
the length of the reattachment zone, L, varies inversely with Re 1 /2 ; A 

therefore, as Re increases, for a given x Q > L, the dependence of the 5 

average heat transfer on Re decreases. 0 

1 

Finally, one can make an approximate comparison of the heat transfer 
at the stagnation point of a hypersonic vehicle with that at a probable 
reattachment zone. This comparison showed that the heat transfer at the 
reattachment zone could be as much as two or more times that at the 
stagnation point when Z/R < 1. 

In an actual separated region, the boundary- layer thickness at the 
separation point will not be zero as assumed in reference 1. For this 
case, the present analysis remains unchanged through equation (22) . The 
heat-transfer equations (25) and (26), however, should be modified with 
the proper values of A and v Q at y = L which may be obtained from an 
analysis similar to that in reference 7- 


Ames Research Center 

Rational Aeronautics and Space Administration 
Moffett Field, Calif., July 10 , 1961 
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APPENDIX 


SOLUTION OF EQUATION (l) BX SEPARATION OF VARIABLES 


In a region defined by 0 < x < L and 0 < y < L the strearawise 
velocity distribution of reference 1 can be closely approximated by 

r ri sinh A(L - x) /. n \ 

v(x, L) = -Vo £ ■ (Al) 

sinh AL 

Following the method outlined earlier the inviscid flow equation (l) 
and Its accompanying boundary conditions can be written for 0 < x < L and 
0 < y < L as 


= A 2 i|r - Av 0 coth AL (A2) 

8x By 

t(x, 0) = 0 (A3) 

t(0, y) = 0 (Ah) 


(L, y) = 0 (A5) 


t(x, L) 


1 ° 

A 


coth XL - 


cosh X(L - x) 
sinh XL 


(A6) 


Condition (12) has been replaced by the requirement that the fluid along 
the wall at x - L be moving parallel to the x axis. That this is the 
case is apparent from the symmetry of equation (l) and condition (4). 

In terms of a new function z, where 


^(x, y) = z(x, y) + ^ coth XL (A7) 

A 

the differential equation and its boundary conditions become 


S Z b Z d.2 

(A8) 

(x, 0) = - — — coth XL 

(A9) 

(0 , y) = - -t— coth XL 
A 

(A10) 
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^ (L, y) - 0 
ox 


z(x, L) 


v 0 cosh A(L - x) 
A sinh XL 


(All) 

(A12) 


This new function can be broken into the sum of three functions such that 


z(x, y) = z x (x, y) + z 2 (x, y) + z 3 (x, y) (A13) 

The superposition principle can be used to break equations (A8) through 
(A12) into three simpler problems. 

For zi: 

V 2 z 1 = A 2 z x 


For z 2 : 


Zi(x, 0) 

zi(Oy y) 


Vo 


coth XX 


= 0 


CiZi 


(l, y) = 0 


Zi(x, L) = 0 


V 2 z 2 - A 2 z 2 
z z {x, 0) = 0 

z 2 (0, y) = - —■ coth AL 

A 

(L, y) = 0 


1 


z 2 (x, L) = 0 


< ITNVO 
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For z 3 : 


V 2 Z3 = A 2 z 3 
z 3 (x, 0) a 0 
23 ( 0 , y) = 0 


Sz 3 

<3x 


(L, 


y) 


0 


z 3 (x, L) 


v 0 cosh A(L - x) 
A sinh AL 


Each of these problems is readily solved by the method of separation of 
variables. Combining their solutions with the aid of equations (AI 3 ) and 
(A 7 ) will give a solution to equation (l). For sufficiently large values 
of AL (AL > 5*3); this solution is equivalent to equation ( 20 ) . That 
this should be so can also be seen by examining equations (Al), (A2), 
and (Ad) which for large AL reduce to equations ( 3 ), (9)> stfih (7)* 
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- Inviscid velocity, pressure, and heat-transfer distributions 

along the vail. 
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